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, Abstract 

Oh ■ We introduce a class of models of semiflexible polymers. The latter are 

■ characterized by a strong rigidity, the correlation length associated to the 

. gradient-gradient correlations, called the persistence length, being of the 

same order as the polymer length. 

We determine the macroscopic scaling limit, from which we deduce bounds 
(N ■ on the free energy of a polymer confined inside a narrow tube. 
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O ' 1. Introduction and results 

■ 

^ ' The purpose of the present work is to introduce and study a family of 

effective models of semiflexible polymers. The latter are polymers endowed 
^ ■ with two characteristic properties: 1) Their thermal fluctuations are governed 

by their bending energy, rather than their tension; in other words, they try 
to minimize curvature rather than length. 2) Their persistence length, which 
can be roughly defined as the correlation length associated to the directional 
correlations between tangent vectors to the polymer, is of a size comparable 
to that of the polymer. Such semiflexible polymers play a crucial role in 
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nature. In particular, the biological function of many biopolymers (such as 
DNA, filamentous actin or microtubules) relies on their semiflexibility, the 
latter providing considerable mechanical rigidity. 

1.1. The model 

The model most often used in the physics literature is the so-called worm- 
like chain. In this model, the polymer is described by a smooth path in 
(higher dimensions are of course possible, but we'll stick to this case in this 
paper), of fixed length 1, which we'll write r(s) with s being the parametriza- 
tion such that the tangent vector t(s) = dr/ds satisfies ||t(s)|| = 1 for all s. 
The energy functional associated to the path is then given by 



where $ is usually taken as $( the parameter k, setting the rigidity 

of the polymer. When the polymer makes only small deviations from the 
horizontal axis, an effective representation of the polymer as the graph of 
a function / : M ^ M becomes possible, associating to a given polymer 
configuration / the energy 



where c is the macroscopic length of the polymer 

The main aim of the present paper is to study discrete approximations to 
such models. Namely, we consider lattice configurations (p in the ensemble 
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equipped with the probability measure 
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where the Hamiltonian T-Cn{^) is defined by 



N 




(2) 
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Here we assume that the parameters and e satisfy Ne — c as — oo, 
where c > denotes the macroscopic length of the polymer, and use the 
standard notation for the lattice difference operators 



(Vv2)fc = - Vk-i , Av9fc (V(Vv2))fe+i = v^fc+i - 2v?fc + v^fc-i • (3) 

Similarly, if the gradient condition on the right end of the polymer becomes 
improtant, we consider the ensemble T^^ (with ^ = {^l, ^r} and d = c?Ar+i) of 
configurations with fixed endpoints and fixed gradients at both extremities, 

^N,d = = iVO = 0, . . . , (fiN+l) ■■ (4) 

fl = ^L, V^Af = + djy+i, (fN+l = d]y+i 

= = (V9o = 0, . . .,'fN+i) : 
equipped with the probability measure 

pLm = pU<fi\<fi^^L)- (5) 



Our aim is to describe the typical behaviour of the trajectories cp E I. 



N,d 

under the measure P^^ from ([5]) with Ne ~ c and e <^ 1. Despite the 
Hamiltonian ([2]) might look unusual, our results in App. |A] show that this 
choice corresponds to the semiflexibility regime, when the persistence length 
and the polymer length are of the same order. In addition, our results in 
Sect. 12] below (see, eg.. Remark 12. II and Remark 12.21) show that the Hamilto- 
nian ([2]) with $(x) ~ as X ~ is essentially the only sensible choice from 
the physical point of view. 

Of course, this model shares the limitations of the macroscopic effective 
model it approximates: It forbids backtracks of the polymer, and the gradi- 
ents of ip have to remain close to zero. For the questions we have in mind, 
these approximations will be harmless. 

Let us nevertheless mention that it is also possible to discretize directly 
the wormlike chain, thus obtaining discrete models of semiflexible polymers 
that are better suited to discuss other properties (for example the effect of an 
external force pulling the polymer, which in general results in a non-trivial 
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macroscopic profile). In particular, there is a natural discrete variant of the 
wormlike chain, in which the polymer is modelled as a chain of hard rods of 
unit-length, with an energy penalizing changes of orientation. This model 
might also be amenable to a mathematical analysis, although this would 
surely generate additional technicalities. 

To our knowledge, the mathematical analysis of models of the form in- 
troduced above is still quite limited. The works closest to ours are 0, [sj, 
in which the effect of an external pinning potential, similar to the problem 
we analyse in Part II is done in the case e = 1, as N ^ oo, both with 
and without a positivity constraint. In particular, it is shown that such 
models display a very different critical behaviour from that for interfaces or 
polymers with tension. Notice however, that setting e = 1 and taking to 
infinity implies that the described polymer is not semiflexible anymore (its 
persistence length being of the order of the lattice spacing, while its length 
becomes infinite). 

Other relevant works deal with the case of membranes, a natural higher- 
dimensional analogue of the one-dimensional polymer considered here. These 
models have also important applications, as they can be used to describe, e.g., 
cell membranes. However, their rigorous analysis is quite involved, and up 
to now only the case of objects of internal dimension at least 4 have been 



successfully studied; see and references therein. 

We finally observe that to simplify our exposition we only discuss discrete 
height models. Similar results can also be obtained by analogous methods for 
continuous height models; of course, there one has to understand the RHS of 

and of similar expressions as the densities w.r.t. the Lebesgue measure. 
The key ingredient of our analysis-the local limit theorem-remains the same. 



We refer the interested reader to the classical monograph [12|, Chap. VII of 



which deals with LLT's both in discrete and continuous setting. 

2. Scaling properties of semiflexible polymers 

2.1. Reduction to the RW case 

The problem above can be reduced to a problem about random walks. 
To do this, consider the process 

a = (V^)., 6= a, (6) 
and observe that, w.r.t. the distribution from ([I])-®, its rescaled increments 
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are i.i.d. random variables. We need to study the large- behaviour of such 
random walks conditioned on the event = —^r- Since in view of (jS]) 

m m— 1 

u = ^i + Yl = ^1 + ^ E ' (7) 

j=2 j=l 

the "gradient" boundary condition in (jl]) reads 

N 

£^Vj = ^N+i -^1 = Vv^jv+i - Vv?i . (8) 
Now, using the relation ([7]), we get 

k k-l 

^k = 'fi + ^^m = k^i + £^{k - (9) 

m=2 j=l 

and rewrite the "right-end" boundary condition from (jl]) as 

N 

£ ^ (AT + 1 - j)r]j = ^N+i -{N+ 1)^1 = ^N+i -{N + l)V<^i . (10) 
i=i 

Clearly, the problem above now reads as the problem of describing con- 
ditional distribution of a RW with i.i.d. steps r], subject to constraints ([H]) 
and (ITOl) . Questions of this type are well understood, see [3] for a recent treat- 
ment of a similar model, so we can simply state the corresponding answers 
and discuss the necessary modifications in the proofs. 

2.2. Functional CLT 

We now turn to the analysis of the fluctuations of the process. In view 
of the applications we have in mind, and the intrinsic limitations of this 
model, we shall restrict attention to "macroscopically bounded" boundary 
conditions (see Theorem 12. II below) . 

For /c = 1, . . . , A^, we consider 

^. = E^., y^ = (11) 



5 



Our basic assumption is (remember that N = c/e) that the first two moments 
of 7] satisfy 

Er/ = , = a%, hm Na% = oo. (12) 

It imphes, for all m = 1, 2, . . . , A^, 



EXm, = EK,^ = , VarXm = ma^ 



N 5 



m(m + 1) 2 \/ m(m + l)(2m + l) 2 

(13) 



2(Ar+l) Q{N + iy 



in particular, the vector Zn ^= {X^, Y^) has zero mean and the covariance 
matrix 

Cov(Z^) = ^ iv(2jv+i) o • (14) 

V 2 "A' 6{N+1) "nJ 

We are going to study the asymptotics of the conditional process 9N(t), 
t G [0, 1], related to the one-point projections 

{Yk I Xn = aN, Yn = Bn) 

with ajsf, bjsf chosen in such a way that the probability of the condition 

P[Xn = ttN, Yisf = 6jv) 

remains positive for all large enough and, for some finite > 0, 

Wn\ + l&jvl ^ 

lim sup ^ — < K . 



More precisely, for t G [0, 1] let 

def 



Nt = [Nt] , (15) 



^ Here and below we use E, Var and Cov to denote the expectation, the variance and 
the covariance of various random variables expressed in terms of the 77- variables. We stress 
that with fixed value of the first gradient ^ and fixed law of the i.i.d. increments 77, the 
probability measure P^(<^) becomes uniquely defined. 
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and define the continuous process 6N{t) via 

/ TTl \ 1 

- = y={Y^\XN = aN,YN = bN) (16) 

at the points t = m/N E [0, 1] with subsequent hnear interpolation for other 
values of t G [0, 1]. Our main result reads as follows: 

Theorem 2.1. Let the independent random variables t] have common dis- 
tribution with variance cr^ satisfying Na"^ ^ oo as N oo. If 

iim = = a , lim = = o , 

then the distribution of 6N{t) converges weakly in C[0, 1] to that of a Gaus- 
sian process 9{t), t G [0, 1], such that for all s, t with < s < t < 1, 

£e{t) =t^{t- l)a + t^(3 - 2t)b , 

e2|'1_+^2 (17) 

Cov(^(s), e{t)) = ^ {2t{l -s)+t~sj. 

In particular, for a = b = we get 

eit) r^M{o,lt%i-tf) , tG[o,i]. 

Remark 2.1. 1. In a sense, the main message of the above result is that 
under sufficiently mild assumptions (i.i.d. increments i] with variance 
cr^ satisfying Naj^ — >• cxd as iV — > oo) the only physically relevant 
potentials $ for the model at hand are convex potentials of the Gaussian 
type, $(x) ~ as a; — > 0. Indeed, for every model satisfying 

Theorem 12.1^ there exists a mesoscopic scale 6 = 6iy such that 
Ns = c/6 ^ oo will still satisfy the condition N^af^ oo. As a result, 
it is possible to discretize our macroscopic polymer so that its behaviour 
on the scale 6 is approximately Gaussian. Consequently, among various 
a priori legitimate choices $(x) ~ \x\", the Gaussian case $(x) ~ 
popular in physics literature, seems most natural for the problems we 
discuss here. 

2. The reader might wish to interpret the limiting process 6{t) as the 
"bridge of the integral of a Brownian bridge" . Indeed, with 

1 
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and X = {Xq, Xi, . . . , X^) satisfying the invariance principle, the scal- 
ing limit of Y becomes the integral of the scaling limit of X, ie., a 
Brownian motion; see the proof of Theorem 12.71 below. Of course, a 
similar interpretation holds for other results in this section. We shall 
leave such observations as an exercise for a motivated reader. 
Also, one might wish to notice that the function m{t) = E6{t) satisfies 

m(0) = m'(0) = and m(l)=6, m'(l)=a, 

which is not surprising since our choice of the exponents 7 and 5 in 
fH3|) guarantees that the limiting process 6{t) shares common gradient 
restrictions with all its discretizations. 

Rewriting 1^ in the form 

ifM,+, = {Nt + l)^i + {N + l)eYN,, 
we observe that the boundary conditions (jlj), 

V'0 = 0, ipi=^L, <^N+l = dN+l, ipN = d-N+l + , (18) 

become 

XN = -{U + U)/e, YN = {dr,+i/{N + l)-^L)/e (19) 
and we can rewrite the theorem above in terms of the "profile process" (p: 
Corollary 2.2. Conditioned on fllSp with the property 

-(a + ^R)/{eaNVN) ^ a , (rf^+i/(iV + 1) - ^L)/{eaNVN) ^ b 

as N ^ 00 (recall that Ne ~ cj the distribution of the process 

{ipN,+i - {Nt + l)Q/{aNVN{N + l)e), te [0, 1] , 

w.r.t. to the probability measure P^rf(-) from converges weakly in C[0, 1] 
to the limiting Gaussian distribution with parameters (fT7|) .- in particular, its 
one-dimensional distributions approach 

J\f(t''{t - l)a + t^{3 - 2t)b, It^l - tf^ 

as N ^ 00. 
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Remark 2.2. In the most popular case considered in the physical literature, 
namely the Gaussian case $(a:) ~ with the Hamiltonian (cf. ([2])) 



the random variables rij = 6~^A(fj have variance a]^ = 0{e~^) = 0{N/c), so 
that by the corollary above the fluctuations of the polymer are of order 

aN^N {N + l)e ~ y/cN = c^^h'^ , 

ie., they live on the macroscopic scale. 

The rest of this section is devoted to the proof of Theorem 12.11 We flrst 
derive convergence of flnite-dimensional distributions of the process ^Ar(t) 
(see Theorem 12.71 below) and then establish tightness of the probability dis- 
tributions of 9^(1) in C[0, 1]. 

We turn now to the proof of Theorem 12. 1[ Let x('u), m G R, denote the 
characteristic function of rj, 

X{u) = Eexp{iuT]} ; 
by the moment assumption above we have: 

x{u/aN) = 1 - y + as M ^ 0. (20) 

Fix an integer A; > and a collection of real numbers tj satisfying 

= to < < ^2 < ■ ■ ■ < 4 < tk+i = 1 . (21) 
Our flrst goal is to prove the central limit theorem for the random vector 

^= 7^ {^N, YNt^ , Ym^ , • • • , YNt. , ^Nt... ) • (22) 

To this end, observe that the corresponding characteristic function reads 

k+l 

Xn{uo,ui,. . . ,Uk+i) = Eexp I — ^-j=(uoXn + ^uiYNt^^Y 
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It is convenient to denote 

fc+i 



4(J) = «o + ElvTT(^*^ + l-^■)^ ^<J<N.. 

1=1 



(23) 



fc+i 

u\x) = uo + ^ui{ti-x)+ , 0<x<l. 
1=1 

Then with u = {uq, ui, . . . , Uk+i^ G ]R*^+^ we rewrite 

N 

x'kiu) = U^{u%{j)/{ar,VN)) (24) 

so that, in view of the asymptotic relation (l20l) and the hmiting assumption 
(T%N — s> oo (as —> oo), we get, uniformly in u from compact sets in M'^"'"^, 

logx^(n) = -— XjK(j)]' + o(l) = -2 + (25) 

By a routine (but straightforward!) induction one deduces the following 
result: 

Lemma 2.3. For s, t with < s < t < 1 denote 

/(i^) = |, 9is,t) = ji3t-s). 



Then for every integer k > the quadratic form 

k+l 



pi k+l 

/ [u\x)fdx= V q 



has the matrix 

^'=m":l (26) 

with the entries 

qoo = 1 , qoi = q/o = f{h) , 1 < / < /c + 1 , 
^h,h = ^hh = gith^th) , l<li<l2<k + l. 



(27) 
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The above lemma together with fl25l) imply the central limit result: 

Theorem 2.4 (Central Limit Theorem). For every fixed k > the dis- 
tribution of the random vector from (1221) converges as N oo to the 
Gaussian distribution with zero mean and the covariance matrix defined 
in (126!) - fl27j) . The convergence of the corresponding characteristic functions 
(123]) is uniform on compact subsets of M'^"'"^. 

Remark 2.3. According to (!22|) . this theorem implies that the fluctuations 
of Xat are of order ai\fVN- Combining this with ([8]) and (fTTj) . we see that 
the end-to-end gradient fluctuations Vv^at+i — Vy^i of the polymer are of 
order eaN^/N. In the natural Gaussian scaling a% = 0{e~^) of Remark 12.21 
this implies that 

Var(Vy^^+i - Vy^i) = ©(eV^iV) = 0(c) . 

In other words, the persistence length and the polymer length in our model 
are of the same order. 

Our next goal is to establish the local version of the above theorem. For 
X = (xo, xi,...,Xk, Xk+i^ e M.''^^ 
let Pq(^) denote the probability density 

p^(^) = (2vr)-(^^^^/VetQr^/2exp{-i(Q-^x,s)} 
of the limiting Gaussian distribution with the characteristic function 

Xq{u) = exp|-|(Qw, 



Theorem 2.5 (Local CLT). Let a sequence of vectors 



X 



(N) _ ( {N) (N) (N) (iV)^T fc+2 



be such that x^'^^ ^ x e M''+2 as N ^ oo and the probability P(Z^ = S^^)) 
be positive for all N large enough. Then as N ^ oo we have 

^fc+2^(3fc+4)/2 p(^fc ^ -iN)^^ ^ pfc ^ 

with the remainder o(l) vanishing asymptotically, as N —>■ oo, uniformly in 
X on compact subsets of M'^"'"^. 
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Proof. The claim of the theorem follows from standard considerations pro- 
vided the off-line property is estabhshed (for a recent exposition, see, e.g., 
0, Thm 4.2]); it thus remains to verify the latter. 

By the assumption on the distribution of rj, we have, for all C > small 
enough 

sup \x{u/aN)\=rxe (0,1), (28) 

(:<HaM\<T 

where T = n/d for lattice distributions of period d > 0. In view of the 
factorization fl2^ . the off-line property shall follow once we show that for 
some small enough C = Cfc > sufficiently many values u^{j) satisfy the 
condition (recall fl23|) ) 

\u%ij)\ > CfcCTjV, 

uniformly in large enough. However, by the very definition fl23|) . the 
sequence u%{i), j = 1, 2, . . . , is a piecewise linear sequence of real numbers 
interpolating the values 

and having increments (recall (fT5l) ) 

^ k+l 

- u%ij + 1) = $^w/%<iv,^ • 

1=1 

By 0, Lemma 4.4] it is enough to show that for ^ > as in (128|) one has 
(recall (USD, ([23D) 

max \u%{NtM > 2((rN , 

l=0,...,k+l 

as then the rest of the proof of Theorem 12.51 would be analogous to that of 
0, Thm 4.2]. 

We prove the remaining condition by verifying the following claim. 
Lemma 2.6. For a fixed collection ti, . . . , tk as in fl2T|) . let 

A = mm(tm+i -tm) > . 

Then for every u G M'^^^ such that 

fc+i 

ll*^L = (EwO'">^^^^ 

1=0 
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and all N large enough at least one of the following inequalities holds: 

\u%{N)\ > 2Car, , \u%{Nt^) - u%{Nt^^J\ > ACa^ 1 = 0,. ..,k. (29) 

Proof. We argue by contradiction and start by assuming that none of the 
inequahties fl29l) holds. Since 

u%{N) =uo + Uk+i/{N + 1) 

and 



l>m 

we deduce that 

I I / 4CaAr I I 8CaN i i . 8C0-7V 

\Uk+l\<—^, \uk\<^-, \ui\<—^ 

and therefore that 



\u 



l^<4CV^+(^^)'(« + l)<(^)'(16* + 5) 



We now deduce convergence of finite-dimensional distributions of the pro- 
cess 9n(-) from f|T6l) : 



Theorem 2.7. Let real sequences a^, he such that 
hm = = a , lim 



and the probability P{Xn = cln, Yn = ^at) be positive for all N large enough. 
Then for every k > 1 the k-dimensional distributions of the process 6n{-) 
converge to those of a Gaussian process 0{-), whose parameters are 

E9(t) = t^{t - l)a + ^2(3 - 2t)h , t G [0, 1] , 

s2ri-t)2 (30) 
Coy{e{s),9{t)) = ' [2t{l-s)+t-s\, 0<s<t<l. 

Proof. As the convergence result follows directly from the local limit the- 
orem, we shall only derive the parameters (!30!) of the limiting process 0{t). 
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To start, fixO<s<t<l and notice that the conditional distribution of 



aNVN 
converges to that of 

{Jis),Jit) \wi = a,J{l)=b), 

where 

nl nv nv 

J{v) = (v — u)^ dwu = {v — u) dwu = / Wudu 



and Wsi s G [0, 1] is the standard Brownian motion (Wiener process). Using 
Lemma 12.31 and the classical property of conditional multivariate Gaussian 
distributions, we deduce that the mean of the limiting process equals 

£{j{t)\w, = a,J{l)=b) = {f{t) g{tA)){^^ Q 

= f{t - l)a + t2(3 _ 2t)h 

and its covariance matrix is 

h{s,s) h{s,t)\ ^ /'g{s,s) g{s,t) 
h{s,t) h{t,t)) \g{s,t) g{t,t) 



f{s) gisA)\ /'I iVV /(.) m 

fit) g{t,l)) [l ij {g{s,l) g{t,l) 
with 

h(s,t) = ^ ^-\2t(l - s) + t - s] , 0<s<t<l. 

6 

It remains to prove tightness of the sequence of probability distributions 
of the processes 9n{-) in the space C[0, 1] of continuous functions on [0, 1]. 
To this end it is sufficient (0, Thm 9.2.2]) to show that for some positive C 
and 7 > 1 the inequality 

^\ON{t) - 9n{s)\^ < C\t - (31) 
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holds uniformly in [s,t\ C [0, 1] and all N large enouglijf] The key to 031 p is 
the following result whose proof shall be postponed till the end of the section. 

Lemma 2.8. Let real sequences a^, &Ar he such that 



lim 



hm = 



and the probabihty P(X7v = aN, Yn = b^) be positive for all N large enough. 
There exists a positive constant Ci such that the inequality 

£{Xl I = ar,, = b^) < C,a%N 

holds uniformly in k = 1,2, N . 

The target condition (131 p is a straightforward corollary of the above 
lemma. Indeed, it follows from ffTTj) . the definition f|T6l) and the lemma that 
for every m = 1, 2, . . . , 



m 



N 



In 



m — 1 
N 



< 



X 



N 



aN,Yj 



N — f^N 



(iV + l)2 

Now, observing that for all s, t with 0<s<t<lwe have 



Nt+l 

ONit)-eN{s)= ^ ajv(j) 

j=Ns+l 



N 



-e 



N 



N 



where ajsiij) = 1 for all j in the sum (with possible exception of the extreme 
values j = Ng + 1 and j = Nt + 1, for which aAr(j) G [0, 1]), the Cauchy 
inequahty gives 



Nt+l 



'J - 1 - 

N , 



^ Actually, our argument shows that here 7 = 2; this is not surprising, as the trajectories 
of the limiting process 6'(-) have continuous derivatives. 
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and thus implies the target estimate (1ST]) : 

E|^^^(t) - < C ^^\l\\y^' <C\t-s\\ 

uniformly inO<s<t<l and all large enough. 

We turn now to the proof of Lemma 12.81 and shall treat separately the 
two cases k"^ > N and k"^ < N. 



Proof (of Lemma EIHD- Case k > VTV. Let k = k^ > VN and k/N 
K, E [0, 1] as — oo. Then for the vector 



U'' ^— [ Tf^k, T^^N, Tf^^N 

the central limit theorem holds. Indeed, by a straightforward computation 
we deduce that the characteristic function of satisfies 

lim log E expjif ^° Xk H Xn H Yn] \ 

= (^^^0 + 2t;ot;iv^ + v\ + 2vqV24^(^ ~ i) ^^^^ 

As the variance of the limiting conditional distribution is 

(l - 4k + - G [0, 1] 

and its mean is boundedjll we deduce that for some C2 > 

1 

uniformly in k under consideration. 

Case k < y/N. Using arguments similar to those in j^, pg. 257], we deduce 
that for all j = 1,2, . . . , 

^(^{VjY Xn = ttN, Yn = &Arj < Cscr^ 

(in fact, as explained in j8j] for large N the LHS is close to Erjj = a^). As a 
result, the Cauchy inequality implies 

^(^{XkY Xn = aN,YN = &iv) < PmaxE(^{r]jy 

<Csa%e <C3CT%N. 
The proof of Lemma 12.81 is finished. 



1 1 1 



Xu 



Xn — aN, Yn — bN] < C; 



Xn — aN, Yj 



AT — WAT 



being a linear combination of the constraints a and b (with K-dependent coefficients) ; 
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2.3. Large deviation regime 

By combining the arguments above with the approach of j^, one can 
also describe the large deviation behaviour of semifiexible polymers. As such 
generalization is straightforward, we only mention some results. 

Let Liq{h) denote the log moment generating function of the step distri- 
bution (recall P), 

def r 6"^*^'')+''^ 

= logEexp{/ir/} =log— jr-:^^(^y^^; (32) 

we shall assume that 

L'^(O) = , L^(0) = 4 G (0, +oo) , (33) 

that Ln{-) is finite in some (in general, e-dependent) neighbourhood of the 
origin, and that Ln{-) behaves properly under rescaling: 

LN(—)^L{h) asN-^oo, (34) 

where L(-) is a strictly convex function in some /i-neighbourhood of the origin. 
E.g., for the Gaussian case $(x) = we obviously have 

L^{h) = —h^ and L{h) = -h\ 

^cKj 2 

2.3.1. Probability of the right-end boundary condition 
Let Xm, Ym be as defined in f fTTl) . 

m 1 
j=l j=l 

and let P denote the probability distribution of the RW with steps rjj] we 
shall assume that the assumptions ( !32l) -( !34l) hold. Then the probability 
of the right-end boundary conditions given the left-end ones (essentially of 
finishing a "droplet" at time with gradient — ^r) is 

P(c^G J^,J<^o = 0,V<^i = ^l) 

= P{vn+i = o-iN + 1), Vv2jv+i = -^R I V^o = 0, Vv^i = ^l) 
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and, in view of the relation 

= {[tN] + 1)^1 + (iV + l)eY[tN] , (35) 
clearly, coincides with the LD-type probability 

p(Xn = -e-\^R + a), Yr, = e-\a - ^l)) • 

Its limiting behaviour is well known (see, e.g., 0, Theorem 4.2]), so we just 
recall the corresponding result: 

For real numbers u and v, denote 

Ln{u,v) = logEexpj— Xat + — Yat) ; 



then, as oo, we have 



N ^Ln{u,v) ^ Lao{u,v) ^= / L(ii + {1 — x)v) dx (36) 

Jo 

with L{-) from (132!) . The optimal tilts u*, v* can be determined from the 
conditions (cf. ^, Eq.(2.26)]) 

L'[u + yv) dy = — ^^ ^ 



a - a ^^^^ 

yL'iu + yv) dy = — , 

c 

where we use the fact that ipi\i/N—>-a as N—>-oo in such a way that Ne —>■ c, 
the macroscopic length of the excursion under consideration. Then the sharp 
LD asymptotics for the probability of interest, up to a factor of (1 + o(l)), is 

exp 



[U*, V* 



n{-^^^^*-^^^*-Loo{u*,v*)'^} , (38) 



where D{u,v) stands for the Hessian of L^o as the function of u, v. Clearly, 
the expression in the exponential is just the convex dual of Loo evaluated 
at the point with coordinates as in the RHS of fl37|) . 

A straightforward computation in the Gaussian case ^{x) = kx'^/2 based 
upon the correspondence ( |T8l) -( |T9|) as well as the moments (fT3l)-( fT^ gives 
the following exact analogue of (!38l) for a = and > 1: 



K I12{N + 1) r (2Ar + 1)^2 -2(Ar + 2)^LeR + (2Ar + 1)^2 ■ 



27rm\ N-1 c{N-1)/{Nk) 
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2.3.2. Mean profile 

To catch the mean profile, fix a real t, < t < 1, and consider the vector 

Since according to the relation ([9]) we have the conditional distribution of 
i.p[tN]+i given (Xat, Y/v) can be directly derived from the local limit theorem 
for the vector Z\f. 

Mimicking jj], we introduce the log moment generating function L]^{u, v, w) 
■t 



of the vector 



def (Hi U XV 

V, w) = log Eexp<^ — Xn H Yn H Yum \ , 

and observe that the conditional mean value of the last component Y^tN] of 
given the value of the first two is, up to a small correction, as 

d 



E{YitN] I Xn = xn,Yn = yisf) ^ aN-^L^N{u\v\w) 



10=0 



with the optimal values u*, v* obtained through an analogue of ([3 

= {xn, Vn) ■ 



^ Tt ! \ ^ Tt t 



(m*,i;*,0) 



Observing that (where for a real x we write x"*" == max(x, 0)) 



L%{u,v,w) = 



N 

^ \(Tn N + 1 aN V N + 1 ) a^) 



we immediately obtain, up to a small correction. 



1 d 



Lir(u\v*w) 



w=0 JO 



{t - L' {u* + {1 - x)v*) dx, 



and thus the (conditional) mean value of the macroscopic polymer at "time" 
t is (recall (^) 



(39) 



t^L + c I {t- x)L' [u* + (1 - x)v*) dx . 
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In particular, in the Gaussian case $( /2, the mean rescaled profile 

dsn]) becomes 

It is instructive to compare the previous results to their analogues for the 
interfaces. Of course, the non-trivial geometry of the mean profile as well 
as anomalous C ^-smoothness of the trajectories (recall the comment to (1311) 
above) are due to the nature of semiflexible interaction and are not present 
for interfaces. 



3. Free energy of a confined polymer 

As an application of the above estimates, we turn now to a problem that 
has often been studied in the physics literature (see, e.g., [13| and references 
therein): Determine the free energy (per unit of macroscopic length) of a 
semiflexible polymer constrained to lie inside a tube of given radius. From 
the mathematical point of view, this is equivalent to studying the logarithmic 
asymptotics of the probability of the event {supo<fc<7v+i \ fk\ < p}, when 
is large enough. 

Using the functional CLT, it would be sufficient to prove the correspond- 
ing claim for the limiting Gaussian process. This so-called small ball problem 



has been studied for the integrated Brownian motion in [10|. We are going 
to give a completely different proof, in the spirit of j^, which is easy and 
more robust, and also holds for positive values of e. 

Theorem 3.1. Let c > be the macroscopic length of the polymer. There 
exist constants po = Po(c) > 0, Ci > 0, C2 < oo and 6 > such that, for all 
P < Po, 

CI/ \ C 

< -- log sup^ l^.l < par^VN | ^0 = V^i = o) < , 

uniformly in e < Sp^/'^c^^^ . 

Remark 3.1. 1. The existence of the limit as e = c/iV can be 
proved using a standard subadditivity argument, see [lo|. An exphcit 
expression for the limit seems to be unknown (although the physicists 
have good numerical estimates). 
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2. A similar result holds for other boundary conditions, as long as (fQ is not 
chosen too close to the boundary of the tube, and ,^1 is small enough. A 
similar remark applies for <fN+i and ^n+i (which were unconstrained 
above). For example, the proof remains unchanged if the boundary 
conditions at both extremities satisfy the same constraints as demanded 
by the event A in the proof. 

3. Although the above expression might look superficially different from 
the one given by the physicists' derivations, they actually coincide. To 
see this, it is best to restrict attention to the case studied in the physics 
literature, in which the Hamiltonian is of the form ^ ^^i(^V^i)^, and 
to write down explicitly the temperature dependence. In that case, 
ajf = N/{Pk.c), where /? = l/{k-QT) is the inverse temperature. To 
match the physicists' procedure, we wish to measure the width of the 
tube in units set by the polymer length. The event we are interested 
in thus becomes 

sup A^~^/^|</?fe| < r/c, 
i<ifc<Ar 

where we have denoted by r = ^J c/ {Pk)p the macroscopic width of the 
tube. We then see that the free energy is given by k-oT (/?/t)~^/^r~^/^, 
which agrees perfectly with the physicists' expression, since is the 
persistence length corresponding to these parameters. 

Proof. Lower bound on the probability. We write 

R = pajvViV and D = [e-^'^a-^'^ R^l'] = [p'/'c'/h''] . 

Let also z/ > be a small number (to be chosen below) and denote by A the 
event that 

• \VkD+i\ < i^R, for all 1 < A; < [N/D]; 

• |6d+i| < i^R/D, for all 1 < A; < [N/D]. 
We then have the lower bound 

Pjv( sup \ipk\ <R\^o = ^i = 0) 

l<k<N 

> Pn( sup \ipk\ < R\(po ^ ^1 ^ 0,A) Pn{A \ (po^^i^ 0). 
i<fe<Ar 
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Let us first find a lower bound for Pn{A\(Po = fi = 0). Conditioning on 
the pairs fkD,'^kD+i, 1 < ^ < W/^] (compatible with the event A), the 
Markov property implies that it is sufficient to consider what happens in a 
single piece {{k — l)D, . . . , kD + 1}. Namely, for |ao| < and \go\ < z/, it is 
enough to prove that 

Pd(Iv5d+i| < i^R, I^D+il < i^R/D I Lfo = aoR,^i = goR/D) 

is bounded away from zero, uniformly in c, p and e. Rewriting this event in 
terms of the random variables X£, and Yd yields 

/ 1 ^\ (i^ + ao + go)R - uq- go)R 



aNy/D + 1 L eD^nVD + 1 sDaN^D + 1 J 

(z/ + 5(o)i? {u-go)R 



^ Xd e 



cnVD + 1 L eDcTiyVD + 1 sDaN^jD + 1 

Since 

-^^-^ = l + o(l), asB^oo, (40) 

the Central Limit Theorem 12.41 implies that the above probability converges, 
as D — > cxD, to 

p{Zi e [-{u + ao + go),iy - ao - go], Zq e [-{i^ + go), - Qo]^ 

>p(ZiG [-3z/,-z/],ZoG [-2z/,0]), 
where {Zq, Zi) is a Gaussian vector with zero mean and covariance matrix 



1 l/2\ 
1/2 1/3 i ■ 



This probability being bounded away from zero, uniformly in e, p and c, we 
conclude that 

Pat^^ I (yjo = V'l = Uj > e ' =e ^ , 

uniformly in e, p, c such that D ^ p^l^c^l^e"^ is sufficiently large. 
Let us now turn to the derivation of a lower bound on 

Pat ( sup Iv^fcl < i?| v^o = V^i = 0,^). 

l<fe<Af 
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For |ao|, ja/j+il < v and \g\\, \gD+\\ < ^, let us introduce the event 

E = B{ao, ttD+i, gi, gn+i) = jv^o = aoR, ipo+i = ao+iR, 

i^=g,R/D,in+i = 9D+iRlDy 

Changing to the X, Y variables yields, 

Pd{ sup 1^,1 >R\B)< P( sup ^== > ^^"^/^ I B). 

Fixing some u < the functional CLT and fj40|) then imply that, for all D 
large enough, the latter probability is bounded above by 

p(sup \e{t)\>l), 

where Q{t) is the Gaussian process characterized by f[T7l) with a = gn+i — gi 
and b = ajj+i — a-o — gi- An application of Fernique's inequality [6|] shows 
that this probability is bounded above uniformly in ao, gi, an+i, go+i in the 
range considered. The Markov property then implies that 

Pjvf sup \ifik\< R\ ^0 = ^1 = 0, a) > e~cp'^i^c^i^ ^ 

^l<k<N ' 

uniformly in e, c, p such that D (?l'^c^l'^e~^ is sufficiently large. This com- 
pletes the proof of the lower bound. 

Upper hound on the probability. As for the lower bound, we partition the 
tube into disjoint pieces of length D = [£:~^/^cr^^''^i?^/^]. We then write 

PAr( sup \ipk\ < -R I v^o = V^i = O) 

l<fc<Af 
[N/D] 

< TT PAr( sup \(pk\ < R I V^O = V'l = 0, sup \ipk\ < R)- 

■^-^ \i-l)D+2<k<iD+l 2<fc<{i-l)D+l ^ 

(If N/D is not an integer, we simply bound the contribution of the last, 
shorter, piece by 1.) We are going to show that each of the remaining terms 
in the product is bounded away from 1, uniformly in e,c,p, provided D is 
large enough. The conclusion will then clearly follow. 



23 



Using once more the Markov property, we see that it suffices to bound 
sup Pd( sup \ipk\<R\^o = aQR,ii = giR/D) 

\ao\<l 2<k<D+l 



91 



= 1- inf Pd{ sup \(pk\ > R\ipo = aoR,^i = giR/D) 

|ao|<l 2<fe<D+l 

<1- inf Pd(|<^d| > i?| <^o = ao^,6 = ^1^/^)- 

|ao|<l 

ai 

We shall now separately deal with the cases \gi\ < M and \gi\ > M, where 
M is some large enough number which will be chosen below. 
First, 

inf PD{\'fD\> R\'fo = aoR,^i = giR/D) 

\ao\<l 
\9i\<M 

can be bounded below by 

\Yd\ ^ {M + 2)R 



'al{D + l) eDaNVD + 1 



and the Central Limit Theorem 12.41 and ( HOl) imply that the latter converges, 
as D —>■ oo, to P(|^i| > M + 2), which is bounded away from by a constant 
depending only on M. 

Second, straightforward computations similar to those done in the proof 
of the CLT yield 

Ed {vd I '■Po = aoR, = QiR/ D) = (ao + gi + o{l))R, 
Varoi^D I = aoR,^i = giR/D) = \{D + l)D{2D + l)ale\ 
We conclude that, when |ao| < 1 and \gi\> M, 

I Ed {ipD I = aoR, 6 = giR/D) I > (M - 2)R. 

Let us write ^/jd = — ^d{^d \ = o-oR,^! = giR/ Chebychev's 
inequality implies that 

inf PD{\fD\<R\<^o = aoR,^i = giR/D) 

\ao\<l 



\9i\>M 



< inf PD{\tlJD\>{M-3)R\^o = aoR,^i = giR/D) 

\ao\<l 
\gi\>M 

1 

- 3(M-3)2' 
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and the latter is smaller than 1/3, provided M > 4. 

A. Heuristic derivation of the model 

We wish to construct a discretized version of the worm-like chain model 
from Sect. Il.li Given a positive e, we associate to the macroscopic polymer 
profile / the discretized polymer configuration 

Vk = e~^f{ke); 

i.e., we discretize the polymer horizontally with step e and vertically with 
step e"' , where the (yet unknown) parameter 7 has to be determined. In 
order to determine 7, we proceed as follows. To each polymer configuration 
= (v'o = 0, . . . , 9?Af+i), with = [c/e], we associate the energy 

N 

HNi^p) sY,He-'Ay,,). (41) 
For a smooth profile / we then have 

so that the macroscopic expression for the energy is recovered, in the limit 

£ 0, 

{^)^eJ2'^ {fije)) ^ / $ {fix)) dx , (42) 
provided the relation 

7 + 5 = 2 (43) 

is verified. The above computation holds for all 7, 5 > satisfying (143|1 . 
Here, we choose 7 = 5 = 1, so that for a sufficiently smooth profile /(■) we 
have Vifk ~ f'i^k), i.e., the macroscopic and microscopic gradients coincide. 
As shown in Sect. [21 for the class of models considered in the present pa- 
per this scaling results in both the vertical fiuctuations and the end-to-end 
gradient-gradient fiuctuations for such polymers being macroscopic. This, in 
particular, implies that the persistence length and the polymer length are of 
the same order. 
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